Modified critical correlations close to modulated and rough surfaces 
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Correlation functions are sensitive to the presence of a boundary. Surface modulations give rise to 
modified near surface correlations, which can be measured by scattering probes. To determine these 
correlations, we develop a perturbative calculation in deformations in height from a fiat surface. The 
results, combined with a renormalization group around four dimensions, are also used to predict 
critical behavior near a self-affinely rough surface. We find that a large enough roughness exponent 
can modify surface critical behavior. 
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Bulk properties, such as magnetization, as well as cor- 
relation functions are modified on approaching a surface. 
In particular, critical behavior near surfaces or defects, 
which is quite different from the bulk, has been exten- 
sively studied theoretically (l]-||], by experiments 
and in simulations ||. Along with this development, the 
method of grazing incidence of x-rays and neutrons 
has become a standard tool of probing critical behavior 
near surfaces and interfaces 

Most theoretical studies have been restricted to flat 
surfaces. However, real surfaces mostly deviate from this 
idealized picture. Possible deviations can be divided into 
two classes: (i) advanced experimental methods, e.g., 
lithographic preparation, allow one to endow surfaces 
with specific, regular geometrical patterns down to the 
nanometer scale, with important applications in technol- 
ogy and material science |J; (ii) surfaces or interfaces 
can be naturally rough, e.g., due to growth, fracture, or 
erosion. Most common are self-affinely rough surfaces, 
in which the root mean square height fluctuations on a 
length scale L grow as L^, where £ < 1 is the so-called 
roughness exponent. Self-affine scaling is predicted by 
many numerical and analytical models of surface growth 
JTo| , and is also observed in a number of experiments [jll| . 

In this work we show that the shape of the surface 
has a distinct influence on the properties of an adjacent 
medium with long-range correlations. This is demon- 
strated for two-point correlation functions near a critical 
point of the medium, for both cases (i) and (ii) outlined 
above. The diffuse scattering of x-rays and neutrons at 
grazing incidence due to such correlations appears in ad- 
dition to what would have been observed if the surface 
was separating two homogeneous media ]T^ ] . The modi- 
fied correlations may thus provide an additional and in- 
direct means of characterizing the surface profile. This 
may be of value when other techniques are not possible, 
as in the case of the interior surface of a glass, or an in- 
ternal crack, whereas scattering from a critical fluid or 
binary alloy coating the surface may be feasible. 



In order to study the effects of the surface shape, we 
develop a perturbative expansion of two-point correlation 
functions in the deformations of the height profile. Ini- 
tially for a Gaussian field, the calculations are carried out 
to second order. Already at the first order, the two-point 
correlation functions track the profile from the substrate, 
with a modulation that decreases with the distance of the 
two points from the surface. This leads to explicit predic- 
tions for the structure factor, as a function of the lateral 
wave vector transfer, for a modulated surface [see Eqs. (Q) 
and (@) below]. For example, for a sine modulation with 
wavelength 2jr/k along one direction, the incident wave 
vector is scattered by k in that direction, with amplitude 
proportional to the modulation in height. 

For self-affinely rough surfaces, second order calcula- 
tions are necessary, as the first order results vanish on av- 
erage. For a massless Gaussian field we find the expected 
result that self-affine roughness leads to subleading cor- 
rections to the decay of two-point correlation functions, 
which at a scale r are smaller by a factor of r -2 ^ 1- ^ 
than the leading contribution coming from a flat surface. 
Typical critical systems, however, are described by a non- 
Gaussian (interacting) field theory. In this case, the cor- 
relations are calculated perturbatively in the strength of 
the interaction, and the results interpreted with the aid 
of the renormalization group (RG) in 4 — e dimensions 
[see Eqs. (|ll|) - ( |l4| ) below]. We find that the subleading 
corrections now fall off with a slower power as compared 
to the Gaussian case and, surprisingly, for a sufficiently 
large £ even dominate, giving rise to novel surface critical 
behavior. 

Fluctuations in the critical system located above and 
bounded by the surface will be described by an n- 
component field $(r) = [$i(r), . . . , 3> n (r)], where r is 
assumed to be d dimensional. For example, n = 1 can 
represent an Ising magnet or binary alloy, while n = 2 
may describe a superfluid. The fluctuations are described 
by the statistical Boltzmann weig ht e-P H , with 
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PH{*} = J d d r |i(V$) 2 + f $ 2 + ^(* 2 ) 2 }, (1) 

where uq is the strength of the interaction, set to zero in 
the Gaussian theory, and to ~ T — T c . The above expres- 
sion must be supplemented by a boundary condition on 
the surface. We choose the Dirichlet boundary condition 
$ = 0, which represents the so-called ordinary surface 
universality class, appropriate to magnets, binary alloys, 
and for 4 He near the normal to superfluid transition point 

In the absence of overhangs and inlets, the surface pro- 
file can be described by a single-valued height function 
h(pz), where x spans a D = d — 1 dimensional base plane 
(see Fig.|lj). The Gaussian correlation (<&i(r)%(r')) = 
Sij G(r; r'), according to Eq. (Q) with u$ = 0, can be cal- 
culated using functional integral methods and is 
given by 

G(nr!) = G b (r;r!) - J d D x J d D y (2) 
x G b (r; x, /i(x)) M(x, y) G b (r'; y, h(y)) . 

Note that we denote d dimensional vectors with under- 
lined letters, and D dimensional vectors with boldface 
letters. Position vectors r are thus decomposed accord- 
ing to r = (rn,z), where rii comprises the D = d — 1 
components parallel to the surface and z is the distance 
from the base plane; G b is the bulk correlation function 
of the Gaussian theory, and the kernel M(x, y) satisfies 

J d D yM(x,y)G b (y,h(y);y',h(y')) = 5 D (x-y'). (3) 

While the above results are generally valid, we focus 
on the behavior of the correlation functions at the bulk 
critical point, i.e., for T — T c , where correlations are 
strongest. 

The solution for G can be expanded in a series Go + 
Gi + G-2 + ■ ■ ■ in powers of h(x). The lowest order re- 
sult, G (r;r') = G b (r\\ , z; rj , z') - G b (r\\ , z; rL -z'), cor- 
responds to a flat surface. The bulk correlation func- 
tion G b (r;r') decays as r _ ( d_2+7? ) for large separations 
r = \r — r'|, where the bulk critical exponent rj is given 
by rj = in the Gaussian theory. In contrast, if both 
points remain close to the surface, Go(r;r') decays as 
r -(d-2+jj||)^ whgrg ^ ; s a sur f ace critical exponent given 
by 77|| = 2 in the Gaussian theory [QJ^]. 

The first order result is given by 

Gi (r; r') = - J d D x A(r „ - x, z) /i(x) A(r|, - x, z') , (4) 

where A(x,z) = J m^ro e' px e~ pz , with p = |p|, has 
the form of a representation of the delta function 5 D (x) , 
i.e.,/ d D xA(x,z) = landlim 2 ^ A(x, z) = 6 D (x). The 
first order result G\ tracks the profile /i(x) of the surface. 
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FIG. 1. Position vectors r = (r\\ , z) and r' = (r|i , z') of the 
two-point correlation function in the critical system located 
above and bounded by a deformed surface. The surface profile 
is described by the height function ft(x) and the vertical dis- 
tances of r and r' from the surface are given by 5 = z — h(r\\) 
and 8' = z' — ft(rj|), respectively. 

For example, for r = \r — r'\ — * oo with z and z' fixed, 
the results for Go and Gi imply the behavior 

G(r;r') - [1 - A{r) - A(r')\ r -(^+v\0 ^ ( 5 ) 

up to terms of order (h/z) 2 and (h/z') 2 . Thus, the lead- 
ing power law is the same as for a flat surface, but the 
amplitude is modulated by the surface deformations in 
the vicinity of ry and r|, by (l6j 

A[L) = %Zl JdP x !*&A{x- nt z). (6) 

Already the results at first order indirectly characterize 
the surface in scattering experiments. To demonstrate 
this, we introduce the Fourier transform of the height 
profile as fc(x) = / (07 e lk x fc(k), with fe(-k) = h(k)*, 
and accordingly the lateral structure factor 5(p, z; p', z') 
corresponding to G(r, r') (where the parallel component 
rii is transformed to the lateral wave vector p). The first 
order result in Eq. (|J) then gives 

S x = - e'P z e- p ' z ' h(p + p') . (7) 

For example, for a sine modulation with wavelength 27r/fc 
along, say, the cc-axis, this implies that the incident wave 
vector component p x is scattered to p' x = p x ± k while 
the other components of p remain unchanged. 

A similar behavior occurs for the contribution to S at 
second order in h, which is given by 

S 2 = - e-P* e-P' z ' J |p - k| h(k) h(p + p' k) . 

(8) 

For the sine modulation, this implies that p x is scattered 
by 2k, 0, —2k. In a scattering experiment with grazing 
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incidence, the length scale perpendicular to the surface 
is set by the depth b the evanescent wave penetrates the 
sample, giving rise to diffuse scattering and thereby prob- 
ing the critical correlations close to the surface || . Since 
this diffuse scattering appears in addition to the contri- 
bution already present away from criticality it can 
in principle be separated out by tuning the temperature 
deviation T — T c . We assume that b is much larger than 
the height of the deformations. In this case, the above 
expansion in the deformations results in an expansion in 
powers of h/b -C 1 for the elastic scattering cross sec- 
tion, which allows one to distinguish the corresponding 
contributions via their intensities. 

The second order results are particularly useful when 
dealing with rough surfaces, where the quench averaged 
first order corrections vanish. In particular, we shall as- 
sume that the rough surface is described by a height func- 
tion /i(x) with /i(x) = 0, and 



[ft(x) - My)] 2 = ^ 2C |x - y| 2 (9) 
d D P 



. e ip-(x-y) p -D+2-2( & -p\ 



where the overbar denotes averaging over different sur- 
face profiles. While at large separations the above corre- 
lations grow as x — y| 2l >, we have also introduced a cutoff 
length A to regulate the behavior of the surface at short 
distances, and an overall amplitude length u. 

A characteristic feature of self-affine roughness is sta- 
tistical translational invariance, since the rhs of Eq. (g) 
depends on the distance x— y| only. This implies that the 
averaged structure factor S is proportional to S D (p+p'), 
and depends on z, z\ and p = |p| only. In order to main- 
tain translational invariance, it is convenient to express 
the results for the correlation functions in terms of the 
local distance 6 = z — h(rn) from the surface rather than 
z (see Fig.jl]). The two-point correlation function must 
now vanish as S or 5' go to zero. In terms of these co- 
ordinates, the leading power law behavior of correlations 
is the same as for a flat surface, but the correspond- 
ing amplitude depends on the roughness and is modified 
by a factor of [1 — a (w / A) 2 ^ 1- '''] as compared to a fiat 
surface, where a > is a number of order unity. The 
subleading correction of order h 2 decays with the sepa- 
ration r = \r — r'\ with an additional factor of r~ 2 ( 1_ f 
compared to the leading term. 

For the interacting field theory, governed by Eq. (|l|) 
with wo ^ 0, standard perturbation theory can be ap- 
plied to get the correlation function 



can then be improved by RG, resulting in power laws in r 
with corresponding surface critical exponents jl].g] . For a 
sclf-affinely rough surface, the expansion of G to second 
order in /i(x) is substituted in the above equation, and 
the quench average is obtained using Eq. (g) . The one- 
loop correction then gives rise to six diagrams of order 
uoh 2 . Similarly as for a flat surface, the new diagrams in 
the quench average at order of u^h? produce logarithmic 
contributions in r, which can again be recast as power 
laws. The final results for the two-point correlation func- 
tion at criticality are summarized below. 

Perpendicular correlations are obtained when r moves 
into the bulk while r' remains close to the surface, which 
implies r = \r — r' | — > oo while 5' is fixed (see Fig. [j]) . In 
this case the correlations decay as 



(*<(r)*i(r')> 
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where the first term corresponds to a flat surface with 
r]± = 1— 5^+§£ + 0(e 2 ). The second term describes the 
effect of self-affine roughness, with an amplitude a de- 
pending on u>, A, and C, and the new universal exponent 

n A- 9 

Jj ± = (2~2() + l-2——e+0(s 2 ). (12) 
n + 8 

Similarly, when both points remain close to the surface, 
i.e., both 8 and 5' are fixed, correlations fall off as 



($i(r)*((r')) 



1 



r rf-2+77|| r d-2+?j|| 



In this case the flat surface is governed by 7711 
e + 0(e 2 ), while self-affine roughness gives 



n+2 
n+S 



n A- 1 

rill = (2-2C) + 2- + 0(s 2 ) 



(13) 

2 - 

(14) 



The corrections due to roughness now decay with a 
slower power as compared to the Gaussian case. In- 
deed, for a sufficiently large roughness exponent £, these 
corrections can even dominate the result for the flat 
surface. The borderline roughness exponent is CI = 
1 — I e + 0(e 2 ) for perpendicular, and a different 



value of £ 
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3 n+2 
2 n+S 



e + 0(e 2 ) for parallel corre- 



lations. This is a surprising result from a naive point 
of view since, due to Q < 1, on larger and larger length 
scales a self-affine rough surface looks more and more like 
a flat surface. Note that this effect becomes only visible 
beyond the Gaussian approximation, which corresponds 

to £ = 0. 



($i{r)$i(r!)) = G(r;r!) - 
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J d d R G(r; R) G(R; R) G(R; r') + O{u 2 ) 



(10) 



For a flat surface, the one-loop addition in u a can be 
regularized and renormalized by minimal subtraction of 
poles in e = 4 — d, leading to logarithmic contributions 
in the separation r = |r — r'\. This perturbative result 



To test the generality of this result, we examined the 
correlations for a d = 2 dimensional XY model below the 
Kosterlitz-Thouless temperature |l7]] . Correlations in the 
spin variables s(r) = e l6 ^ decay with power laws in this 
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case. We found that the surface correlations fall off with 
the simple relative factor of r~ 2i ~ 1 ~ ( '\ We attribute this 
to the Gaussian nature of the fluctuations in the phase 
angle 6(r), which are retained in the asymptotics of cor- 
relations for s(r). 

Thermodynamic quantities and correlation functions 
can be obtained from derivatives of the free energy with 
respect to magnetic fields. To discuss surface behavior, 
we introduce distinct fields hb and h s in the bulk and 
close to the surface, respectively. Assuming that our un- 
derlying assumption of the validity of an expansion in 
h(x) holds, the results for the two-point correlation func- 
tion are consistent with the following form for the scaling 
of the leading singular part of the surface free energy per 
projected area, 



(sing) 



'/>;/».;•" ./'.r ) (15) 



where £ ~ \T — T c \~ u is the correlation length that di- 
verges at the critical point. The first term in square 
brackets corresponds to a flat surface with y b and y s de- 
scribing the relevance of bulk and surface fields, respec- 
tively |QD . The second term gives the effect of surface 
roughness, with £ _2 ( 1_ f) reflecting the average increase 
in area. However, to regain the results in Eqs. dill) - (111) 
we have to use a value of y s = 1 
is different from w, = 1 2_ e 
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e + 0(e 2 ) which 
To motivate and 

justify this difference, we resort to an analogy in which 
the rough surface is replaced with a collection of edges 
with a (possibly scale-dependent) distribution of opening 
angles. Already for a single edge, describing correlations 
requires a distinct and angle-dependent value of y e for 
the magnetic field close to the edge Jl8|,|l9| . Similarly, re- 
sults obtained recently for correlations in the vicinity of 
a fractal surface with fractal dimension df |20|]2lfl cannot 
be obtained using the value of y s for a flat surface [with 
£- d f replacing in Eq. ( [ToT) and omitting the second 

term in square brackets]. Thus y s can be regarded as 
inherently related to self-affine geometry. Interestingly, 
however, it does not depend on the roughness exponent 
C, at least to order e. 

Future extensions of this research could focus on other 
surface universality classes, the possibility of multifractal 
correlations, Casimir-type effects, and in particular exact 
results in two dimensions. However, our results derived 
here already allow for meaningful tests by simulations 
and experiments. 
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